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Abstract Using direct-write atom lithography, Fe na- 
nolines are deposited with a pitch of 186 nm, a full width 
at half maximum (FWHM) of 50 nm, and a height of up 
to 6 nm. These values are achieved by relying on geo- 
metrical coUimation of the atomic beam, thus without 
using laser coUimation techniques. This opens the way 
for applying direct-write atom lithography to a wide va- 
riety of elements. 

PACS 81.16.Ta: 37.10.Vz; 3.75.Be 



1 Introduction 

In direct-write atom lithography, also known as laser- 
focused deposition, atoms are focused into a periodic 
pattern by a standing light wave. This technique has 
been apphed to Na iJj, Cr 013], Al jj, Yb ^Sj, and Fe [i 
[7]. Just as for the case of focusing light, a well coUimated 
atomic beam is crucial for atom lithography. An atomic 
beam arriving at the lens with a local angular spread 
a (see Fig. [T|) will result in broadening of the deposited 
features. To obtain a well coUimated atomic beam, laser 
cooling techniques are generally used. However, effective 
laser cooling requires a closed optical transition. For nu- 
merous elements a closed transition is not accessible with 
present lasers, or hyperfine splitting is present, and the 
use of several repumping lasers becomes a necessity. For 
Fe, for instance, an accessible closed transition from the 
ground state does not exist, resulting in the loss of fo- 
cusable atoms when applying laser cooling [8]. 

It is because of these drawbacks that we decided to 
laser focus Fe without the use of any laser coUimation 
techniques, simply relying for this part on geometrical 
coUimation. We will show that in this way periodic grids 
of Fe nanolines with a FWHM of 50 nm, a height of up 
to 6 nm, and a period of 186 nm can still be deposited. 
Thus, we prove that atoms can be focused without us- 
ing laser coUimation. This makes the technique suitable 
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Fig. 1 Geometrical coUimation. At each spot on the sub- 
strate, atoms arrive with an angular spread a, defined by the 
nozzle size and the distance between source and sample. 



for direct-write atom lithography of any species with 
transitions of sufficient strength that are not necessarily 
closed. Especially for atom lithography of technologically 
relevant elements such as Ga [9 , In [TOlITT] . and Si [12], 
for which laser coUimation imposes large experimental 
difficulties, this is a step forward. Also, laser focusing of 
molecules may be feasible this way. 

For an earlier review of both direct-write and resist- 
based atom lithography, see Ref. [13]. 

As the experiments are time-consuming, it is impor- 
tant to choose the experimental parameters carefully a 
priori. This is facilitated by the availability of an efficient 
and realistic simulation procedure for the atomic trajec- 
tories. Also, a comparison with simulations is essential 
for the careful analysis of the experimental results. We 
use a fast semiclassical Monte-Carlo procedure based on 
the dressed state model [14 that includes the effects of 
diffusion induced by spontaneous emission as well as of 
non-adiabatic transitions between dressed states. 

In Sect.[2]our procedure for simulating standing- wave 
laser focusing is described in more detail. In Sect. [3] the 
experimental setup is described. The experimental re- 
sults of atom lithography without laser coUimation are 
compared with the outcome of the simulations in Sect. 
|4] followed by the conclusions in Sect. [5) 



2 Simulations 

Our simulation is based on the dressed state model in- 
troduced by Dalibard and Cohen- Tannoudji \T^. In this 
model one considers the eigenstates of the complete atom 
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plus laser light field system. The energies of these states 
are given by: 



En - = (n + l)hujL - 



hA nn 

~2 2~' 



(1) 



with n the number of photons in the light field, ujl the 
frequency of the light field, A — ujl — the detuning 
of the light field from the atomic resonance loq . The fre- 
quency separation between the dressed states Q is given 
by: 



n 



A\ 



(2) 



with uji^ = r^yi /2Is the Rabi-freqency, where F repre- 
sents the natural linewidth, / the light intensity, and 7s 
the saturation intensity. The dressed states are a linear 
superposition of the ground and excited states: 



|n,+)= cos(0)|e,n) +sm{9)\g,n + I), 
\n, — ) = — sin(0)|e, n) +cos{6)\g, n + 1). 



(3) 



Here 6 is defined by: 



causing an instant reversal of the force on the atom. The 
relevant transition rates are given by: 



r+_ = Tcos"* 
r. = Tsin'^ 



(6) 
(7) 



with 6 defined in Eq. |4] The probability of a dressed 
state transition due to spontaneous emission is given by 

tF^ or tF y, depending on the initial dressed state, 

with r the time interval. 



< 



S; 



AAAA 




■ « 







0.505 
0.500 
0.495 
-0.495 
-0.500 
-0.505 



,0 0.5 1.0 1.5 
X/A 



2.0 0.0 0.5 1.0 1.5 2.0 
x/A 



Fig. 2 Left: Dressed state potentials for aJij/|Zi| — 10. 
Dashed lines indicate zero-intensity levels. Note the differ- 
ence in shape between the potential minima for |-|-) and |— ) 
states. Right: Same, for aj_R/|Zi| — 0.1. In this case both po- 
tentials are equivalent. 



cos(26l) = -A/Q. (4) 

In our simulation, atoms move as classical point par- 
ticles in a potential field given by the dressed-state en- 
ergy shift zb^^^^Y^i, where the sign indicates the dressed 
state the system is currently in. The light mask is de- 
scribed as a one-dimensional Gaussian standing wave 
along the a;-axis with a waist w and maximum on-axis 
intensity Iq- 

/(r) = lo sin^kx) exp(-2^). (5) 

Here k = 2tt/X is the wavenumber of the light and the 
propagation direction of the atoms is assumed to be 
along the z-axis. As the atoms' kinetic energy is much 
larger than the potential height of the light mask we can 
neglect its effect on the longitudinal motion, so that the 
atom's motion is described by a one-dimensional New- 
tonian equation of motion in the x-direction, which is 
solved numerically. 

Outside the standing wave the intensity I of the light 
is zero and the ground state is |-|-) for Zi > and |— ) 
for A < 0. Unless the atoms undergo a spontaneous or 
non-adiabatic transition, they will stay in their initial 
state. 

The effects of spontaneous emission, and hence dif- 
fusion in the atom's motion, are included in the model 
by evaluating the probability that an atom has under- 
gone spontaneous emission at regular time intervals. By 
far the largest contribution to the diffusion is caused 
by emissions that change the sign of the dressed state, 



Non-adiabatic transitions can occur when the atom 
moves too fast, or in other words, the eigenstates (dressed 
states) change too fast for the internal state of the atom 
to follow. The atom is then transferred from one dressed 
state to the other. This is most likely to occur when 
atoms cross a node of the standing wave. Non-adiabatic 
transfer is incorporated in the simulation procedure, but 
does not contribute significantly for present experimen- 
tal conditions. 

Before we show the results of these semiclassical sim- 
ulations, we discuss the potential that is used. The left 
part of Fig. |2] displays the potential for the case of a 
strong light field, where Lt;/j/|Z\| = 10. The case of a 
weak light field with w^j/jZil = 0.1 is displayed to the 
right. Both potentials are very different in shape. 

In the strong field limit, the potential minima of the 
|-|-) state differ radically in shape from those of the |— ) 
state. The minima of the |— ) state are smooth, broad and 
parabolic over a large range around the minima. Using 
such a potential is a way to avoid non-parabolic aberra- 
tions in the atom focusing process. On the other hand, 
the minima of the |-|-) state look like a V-shaped poten- 
tial. The oscillation period of the atoms in the potential 
then depends on the starting point of the oscillation, 
leading to strong aberration in the focussing. The opti- 
cal equivalent of such a lens is a conical lens or axicon. 
The aberration leads to a less sharply defined focus with 
a long focal depth, or a focal line. 

For positive detuning, the atoms start out in the |+) 
state and are subject to the V-shaped potential. Pro- 
vided that they do not switch between dressed states due 
to spontaneous emission, this then creates focal lines - 
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aligned with the nodes - rather than focal points. Hence, 
the formation of structures should be fairly insensitive 
to the experimental parameters. For negative detuning, 
the situation is reversed and the atoms are focussed by 
the smooth potential, in principle creating sharper fo- 
cal points - aligned with the antinodes - because of the 
reduced non-parabolic aberration. However, these focal 
points are subject to severe abberation caused by the 
spread in longitudinal velocity, the atom-optical equiva- 
lent of chromatic aberration. This negates the advantage 
of the reduced non-parabolic aberration. 

The potential for a weak standing wave is sine-like. 
The difference in shape between the |-|-) and |— ) state 
minima has vanished, and focusing is expected to pro- 
ceed in a similar fashion for both states. 

The experiments and simulations described here have 
all been performed in the strong field case, with typical 
maximum values of u;ij/|Z\| ~ 5, with both positive and 
negative detuning. 

For the experiment, the — s-^Fs atomic transition 
of iron is used. The wavelength is 372 nm. The isotope 
^^Fe has a natural abundance of 91.8%, of which at a 
typical evaporation temperature of 2000 K rv50% is still 
in the ground state. Hyperfine structure is absent 
in ^'^Fe. For this transition, T = 1.62 x lO'^s^^ and the 
saturation intensity for the Mg = 4 ^ = 5 mag- 
netic subtransition is Is — 65W/m^. The light mask is 
a one dimensional standing light wave produced by an 
elliptical Gaussian beam with waist radius Wx along the 
atomic beam axis and Wy perpendicular to it. 

For the calculation, we start with 10,000 atoms. Each 
atom is initially assumed to be in a random magnetic 
substate. The saturation intensity is then calculated us- 
ing the appropriate Clebsch-Gordan coefficient for the 
substate and laser polarization. Spatially, the atoms are 
homogeneously distributed over a single wavelength in 
the x-direction. 

The longitudinal velocity distribution of the atomic 
beam is that of an effusive beam, i.e., a Maxwell-Boltz- 

mann distribution with average velocity {v) = \J~^^^- 
The transverse velocity distribution of the atoms is de- 
termined by the longitudinal velocity and the geometri- 
cal coUimation. The beam divergence, emerging from a 
round nozzle with diameter Z? at a distance L from the 
standing wave, is characterized in good approximation 
by a Gaussian angular distribution with a root-mean- 
square (RMS) width of jj-. This leads to a Gaussian 
transverse velocity distribution with an RMS spread of 

The equation of motion of every atom is integrated 
over a set time interval, before we check for spontaneous 
and non-adiabatic transitions. The effect of both transi- 
tions is to change from one dressed state to the other. 
The possibility that atoms, which undergo a spontaneous 
emission, decay to a different magnetic substate, is ne- 
glected. 
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Fig. 3 Flux distribution of an atomic beam passing from 
below through a strong light mask ujii/\A\ = 10. Left: A < 
generates focal points, i.e., the size of the focus is approx- 
imately one beam waist in the z-direction. Right: A > 
results in focal lines, i.e., the focus is spread over more than 
two beam waists. This difference in flux distribution is due 
to the difference in potential as shown on the left side of Fig. 

El 

The integration starts when the atom's longitudinal 
position (z) is three times the waist radius of the Gaus- 
sian laser beam before the center of the laser beam, and 
ends at the same distance after the center. We make a 
histogram of the atomic flux distribution at set longitu- 
dinal positions. 

To demonstrate the difference in focusing between 
the negative and positive detuning, we first simulate fo- 
cusing of a perfectly coUimated atom beam, setting the 
RMS spread in transverse velocity to zero. A round laser 
beam with Wy ^ = 50 /im, a power of 50 mW and lin- 
ear polarization is taken. The detuning of the light with 
respect to the atomic transition is set at A = ±150 
(27r) MHz (±58 F). These parameters lead to a maxi- 
mum value of cl;/j/|Z\| « 5.9 The results are shown in 
Fig. |3] Interaction with a red-detuned standing wave 
generates focal points, as shown on the left. The blue- 
detuned standing wave gives rise to focal lines, as shown 
on the right. For blue detuning, the atoms experience on 
average 0.45 dressed state-changing spontaneous emis- 
sions occur. For red detuning, the average number of 
emissions is 0.7, reflecting the fact that atoms are at- 
tracted to intensity maxima in this case, where the ab- 
sorption rate is higher. 

The strength of the lenses varies per atom. This is 
due primarily to the distribution over the magnetic sub- 
states of the Fe atoms; secondarily, the longitudinal ve- 
locity spread of the Fe atoms also contributes. The rela- 
tively long z-range over which atom focusing occurs for 
the blue detuned case stands out. 

3 Experimental setup 

An Fe atomic beam is produced using an AI2O3 cru- 
cible with a nozzle diameter of 1 mm heated to a tem- 
perature of around 2000 K by a carbon spiral heater 
[151 , resulting in a typical Fe density in the source of 
npe = 4 X 10^° m~'^. At this temperature the Fe-beam 
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intensity is Ipe = 2.5 x 10^® s ^sr ^, and the average 
longitudinal velocity {v) — 870 m/s. 

The 372 nm light needed to excite the — >^F5 
transition is produced by a titanium-sapphire laser, fre- 
quency doubled with an LBO-crystal in a ring cavity. 
This laser system is locked to the transition using a hol- 
low cathode discharge cell, described in [16^. An Acousto- 
Optical Modulator (AOM) is used to introduce the de- 
tuning oi A — ±150(27r) MHz between the light and the 
atomic resonance. 



Atomic beam 




Fig. 4 Left: sample holder. Right: cross section of sample 
holder. The Si sample is clamped to a ceramic isolator that 
is screwed to the steel frame. The mirror is pressed to the 
frame by a spring. 



The Fe atoms are deposited on a native oxide Si [100] 
substrate, 650 mm downstream from the nozzle. With 
the 1 mm nozzle hole diameter, the divergence of the 
Fe beam therefore amounts to a^MS — 0.4mrad RMS 
and the spread in transverse velocity to 0.35 m/s RMS. 
The divergence is thus considerably larger than the best 
value obtained in our earlier experiments with transverse 
laser cooling of Fe [8] [urms = 0.17m/s). However, the 
divergence does not yet limit the obtained width of the 
deposited nanolines, which is equal to the width obtained 
with the lasercoUimated beam [6 . 

The Si samples can be 3 to 8 mm wide and 2 to 10 mm 
long. The Si sample is clamped on a compact sample 
holder (Fig. [4]), on which also an 8 x 8 x 3 mm mirror 
is mounted. These sample holders can easily be removed 
from the deposition chamber and stored in a load lock 
by a magnetic translator so that the alignment of sample 
and mirror, which have to be perpendicular to each other 
within 1 mrad, can be done ex vacuo. The pressure in 
the deposition chamber is typically 1 x 10~® mbar. 

To create the standing light wave, the laser beam is 
focused to a waist of 50 /im parallel to the atomic beam 
and 80 /im perpendicular to the atomic beam, located 
at the surface of the sample holder mirror. The light is 
linearly polarized. 

The substrate is positioned close to the center of the 
laser beam. Experiments have been performed with the 
substrate at the center of the waist, thus cutting the 
laser beam in half, as well as with the substrate posi- 
tioned at the rear end of the laser beam waist, such that 



just 10% of the power in the laser beam is cut off by the 
substrate. Both cases are schematically drawn in Fig.|5] 
According to the simulations, the 50% cut-off configu- 
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Fig. 5 The cut-off of the laser beam by the substrate. Left: 
The laser beam is cut in half by the substrate (50% cut-off). 
Right: Only 10% of the laser beam is cut off by the substrate 
(10% cut-off). 

ration should lead to somewhat narrower focussed lines 
and results that are less sensitive to variations in laser 
power. However, diffraction effects, which are difficult to 
take into account in the calculations, are more serious 
than in the 10% cut-off case. 

To prevent Fe from oxidizing, a Ag capping layer of 
approximately 5 nm is deposited using an effusive source 
operated at a typical temperature of 1140 K, resulting 
in a deposition rate of 0.15 nm/min. 
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Fig. 6 Left: SEM-scan of the SiN mask. The pitch of the 
lines is 744 nm. Right: AFM-scan of the deposited Fe-lines 
through the mask. Deposition time is two hours. 

Although the height of the deposited structures can 
be easily measured by an AFM microscope, in order to 
determine the average thickness of the deposited Fe layer 
as well (and hence the thickness of the background layer) 
we need to calibrate the overall Fe deposition rate. To 
achieve this, we first produce structures by depositing 
Fe through a mechanical mask. The mask is made by e- 
beam lithography of a SiN membrane. Lines have been 
etched in the membrane with a width of about 150 nm 
and a period of 744 nm over an area of 250 x 250 /im. 
The mask is pressed onto the substrate with a metal foil 
of 100 /im as spacer between the bottom of the mask 
and the front of the substrate. Deposition through such 
a mask results in background-free nanolines, and thus 
gives a direct measurement of the total deposited layer 
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thickness. A SEM-scan of the mask and an AFM-scan of 
the Fe-hnes grown through the mask are shown in Fig. 
m The FWHM width of the structures is about 150 nm. 
From the height of the deposited structures the depo- 
sition rate is measured. The measured deposition rate 
varies between 3 and 7 nm/hour, depending on the exact 
temperature and fiUing of the crucible. In the evaluation 
of the deposited structures, the calibration of the depo- 
sition rate has been performed under conditions as close 
to the specific deposition experiment as possible. 



4 Results 

We have deposited nanoline arrays of Fe with both neg- 
ative and positive detuning of the light and with sub- 
strates positioned both at 50% cut-off and at 10% cut- 
off. The dependence on the light intensity of the height 
and width of the nanolines has been investigated, as well 
as the ratio of line height to the average layer thickness. 
Since we have an AOM with a fixed RF-frequency of 
150(27r) MHz, our measurements are limited to blue and 
red detuned focusing at this frequency shift from reso- 
nance. 

On each deposited sample, the height and width of 
the nanolines has been measured as a function of position 
transverse to the laser beam axis (i.e., along the y-axis). 
Fig. [7] depicts the geometry. As the laser beam has a 
Gaussian profile, each y-position is characterized by a 
different light intensity. 
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Fig. 7 Deposition assembly viewed from the direction of the 
incoming atomic beam. Nanolines are deposited on the sub- 
strate where the atomic beam intersects the standing light 
field (round spot). The height and width of the nanolines are 
investigated as a function of position along the nanolines. 
Since the laser beam is Gaussian, each position corresponds 
to a specific light intensity. 



In Fig. |8] a 4 X 1.5 fim AFM scan of a typical sample 
after deposition is shown. In 2 hours deposition time 
we grow structures up to heights of 6 nm. Fig. |9] shows 
a cross section through the scan over a range of 1 /im. 
The distance between consecutive lines (the pitch of the 
modulation) is equal to A/2=:186 nm. 

Fig. [To] shows a typical profile of the height of the 
nanolines as a function of transverse position. The in- 
tensity profile of the laser beam is shown as well. For 
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Fig. 8 Example of a 4 xl.5 ^m AFM-scan of the nanolines. 
The height (grey)scale ranges over 6 nm. 
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Fig. 9 1 nm line-scan, averaged over 40 nm in the Y- 
direction, of the same sample shown in Fig.jS] The deposition 
time is 2 hours. The height of the structures in this image 
is 6 nm and their FWHM width approximately 60 nm. The 
laser power is P—50 mW, 10% of the laser beam is cut-off by 
the Si substrate. 



high intensities, the height of the nanolines is largely in- 
dependent on the light intensity both in experiment and 
simulation: the height profile does not follow the Gaus- 
sian shaped laser intensity profile, but it has a broad 
fiat top. The height of the potential shown in Fig. [2] 
follows, in the high intensity limit, from the square-root 
of the intensity. Increasing the intensity will move the 
focus over the z — point towards negative z- values. 
Since the lens, especially in the blue detuned case, has a 
large focal depth, the atoms remain focused at the z — 
position, even for higher intensities, resulting in the fiat 
top of the height profile. 

The results of the simulations are shown as well in 
Fig. 



10 



Both curves are normalized to the same height. 
The agreement between the experimental data and the 
results of the simulations is good when we look at the 
shape of the profile. At low intensities, the measurements 
deviate somewhat from the simulations. This may be 
caused by the non-Gaussian tails in the spatial profile of 
the laser beam used in the experiment. 
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Fig. 10 Typical iieigiit profile as a function of position along 
the lines for a blue detuned 50% cut-ofT laser beam. Both sim- 
ulation data and experimental data are normalized to their 
maximal height. The squares are the measurements, the tri- 
angles are the simulation. The patterned profile reflects the 
dependence of the standing wave intensity on position. The 
laser power is P=50 mW. 



Experimental data and simulation results on the width 
of the structures as function of the position on the sub- 
strate are shown in Fig. |11| For high intensities, blue de- 
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Fig. 11 FWHM of the lines as a function of position along 
the lines for blue and red detuning. The squares are the mea- 
surements. Simulations are indicated by the triangles. In the 
lower right graph, the laser intensity profile as used in the 
simulations is shown. 



tuning results in slightly narrower lines compared to red 
detuning in the experiment. This difference is more pro- 
nounced in the simulations, and results from different fo- 
cusing potentials (Fig. [2]). The match between measured 
and simulated widths is better in the 10% cut-off config- 
uration than in the 50% cut-off configuration. Diffraction 
of the laser beam is not included in the simulation. This 
may result in an underestimate of the nanoline widths 
in the 50% cut-off simulation. 

Determining the ratio of the height of the nanolines 
to the average thickness of the Fe layer requires the cal- 



ibration of the overall deposition rate as discussed in 
Section |3] Although this calibration is not very accu- 
rate, the result is that in our experiment the ratio for 
all experimental conditions is approximately 0.5 for the 
lines at maximum laser intensity (the center of the height 
profile as shown in Fig. 



10) 



Table 1 Simulation results of the ratio of the height of the 
lines to the average layer thickness. 





red 


blue 


50% cut-off 


1.25 


1.25 


10% cut-off 


0.60 


0.70 



To obtain the same ratio from the simulations we 
assume that 50% of the atoms are in the ^D4 ground 
state, which is the thermal occupation of the ground 
state at a source temperature of 2000 K. The other 50% 
of the ^^Fe atoms, as well as the other isotopes, are not 
affected by the focusing light field and only contribute 
to the background layer. In Table [T] the ratio between 
simulated line heights and the average Fe layer thickness 
at maximum laser intensity is listed. For the 10% cut-off 
setting the simulated ratio is close to the experimental 
value. For the 50% cut-off case, the experimental ratio is 
much lower than the simulated ratio. This discrepancy 
between measurement and simulation for the 50% cut-off 
case is also present in the widths of the lines: The width 
of the simulated nanolines is smaller than the width of 
the deposited nanolines. This means that the quality of 
the deposited nanolines in the 50% cut-off case is lower 
than what can ideally be achieved. This may be due to 
the diffraction of the laser beam in the 50% cut-off case, 
which is not included in the simulations. 

An analysis of the magnetic properties of the de- 
posited nanolines will be published elsewhere jl7j . 



5 Conclusions 

By direct-write atom lithography Fe-nanolines are de- 
posited with a pitch of 186 nm and a FWHM width of 
50 nm, without the use of laser coUimation techniques. 
Sufficient coUimation is obtained by strong geometrical 
coUimation, produced by the relatively large distance 
(650 mm) between Fe-source and sample compared to 
the 1 mm source nozzle. In this way the divergence of 
the atoms arriving at the sample is reduced to 0.4 mrad 
RMS. 

Experiments and simulations are performed with the 
focusing standing light wave cut in half by the substrate 
and with a 10% cut-off by the substrate. Experiments 
and simulations are in good agreement with respect to 
the height of the nanolines. The experimental width of 
the nanolines corresponds well with the simulations for 
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the 10% cut-ofF case. For the 50% cut-ofF case the mea- 
sured width of the hnes is shghtly larger than the simu- 
lated value. Also, the measured height of the lines versus 
average deposited layer thickness ratio is lower than in 
the simulations in that case. We attribute this discrep- 
ancy to diffraction of the light beam by the substrate 
and small misalignments of the standing light wave. 

The fact that direct-write atom lithography can be 
applied without the use of laser cooling techniques opens 
the way to many new applications. Except for a small 
number of elements with a strong closed-level transition 
in an easily accessible wavelength region, laser cooling 
can add a serious and sometimes insurmountable com- 
plication to atom lithography. Past efforts to apply the 
technique to atoms that are interesting from a technolog- 
ical viewpoint (e.g., gallium or indium in view of III-V 
semiconductor applications) have concentrated on the 
laser coUimation as a first step. These efforts have suc- 
ceeded in the challenging goal of achieving laser coUi- 
mation |51ITT|. but have not yet resulted in the produc- 
tion of nanostructures. The presented results show that, 
by using an high-flux beam source and relying on geo- 
metrical coUimation, laser coUimation can be bypassed 
completely. 
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